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ABSTRACT 

Up to now, most stellar-mass black holes were discovered in X-ray emitting 
binaries, in which the black holes are formed through a common-envelope evolu- 
tion. Here we give evidence for the presence of a massive black hole candidate as 
a tertiary companion in the massive eclipsing binary V Puppis. We found that 
the orbital period of this short-period binary (P=1.45 days) shows a periodic 
variation while it undergoes a long-term increase. The cyclic period oscillation 
can be interpreted by the light-travel time effect via the presence of a third body 
with a mass no less than 10.4 solar mass. However, no spectral lines of the third 
body were discovered indicating that it is a massive black hole candidate. The 
black hole candidate may correspond to the weak X-ray source close to V Puppis 
discovered by Uhuru, Copernicus, and ROSAT satellites produced by accreting 
materials from the massive binary via stellar wind. The circumstellar matter 
with many heavy elements around this binary may be formed by the supernova 
explosion of the progenitor of the massive black hole. All of the observations 
suggest that a massive black hole is orbiting the massive close binary V Pup- 
pis with a period of 5.47 years. Meanwhile, we found the central close binary 
is undergoing slow mass transfer from the secondary to the primary star on a 
nuclear time scale of the secondary component, revealing that the system has 
passed through a rapid mass-transfer stage. 

Subject headings: Stars: binaries : close - Stars: binaries : eclipsing — Stars: 
individuals (V Puppis) — Stars: evolution — Stars: black hole candidate 
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1. Introduction 

With an orbital period of 1.4545 days and two B-type components stars (spectral types 
Bl and B3), V Puppis (HR3129 and HD 65818) is one of a few massive semi-detached 
eclipsing binaries where the less massive components are filling their Roche lobes (e.g., 
Andersen et al. (1983); Bell et al. (1987a, b); Terrell et al. (2005)). According to the 
theory of binary evolution (e.g., Sybesma 1986), theses systems were formed through a 
case A mass-transfer (mass transfer is occurring when hydrogen is burning in the core). 
Since the discovery of V Puppis more than one century ago, it was extensively investigated 
photometrically and spectroscopically (e.g. Andersen et al. (1983) and references therein). 
A weak X-ray source in the field of V Puppis (e.g., Giacconi et al. (1974); Groote et al. 
1978); Bahcall et al. (1975)) and a significant amount of circumstellar matter around the 
binary (e.g., York et al. 1976 and Koch et al. 1981) make it a very interesting binary to 
study. 

Because of its stable light curve and the deep and sharp eclipsing minima (e.g., Andersen 
et al. (1983)), times of light minimum of V Puppis can be determined precisely. Thus, the 
variations of the orbital period can in principle be derived with high precision by analyzing 
the observed-calculated (O-C) diagram. In the present letter, orbital period changes of V 
Puppis are investigated. Then, based on the period variations, the presence of a massive 
black hole companion and the evolutionary state of the system are discussed. 

2. Analysis of the orbital period changes 

Epoches and orbital periods of V Puppis have been given by several authors (e.g., 
Andersen et al. 1983). Kreiner & Ziolkowski (1978) did not detect any period change for 
the binary, while Andersen et al. (1983) pointed out that the period of V Puppis is variable. 
However, since the early eclipse times were mainly observed photographically or visually, 
conclusions on the period change are not reliable. For the present analysis, we only use the 
photoelectric (Pe) and charge-coupled device (CCD) observations. Most of the data were 
collected by Kreiner et al. (2001). Two new CCD times of light minimum were obtained 
on November 2, 2007 and on March 24, 2008 with the Virpi S. Niemela 0.8-m telescope at 
La Plata Observatory (see Table 1). During the observation, the V filter was used. The 
(O — C)i values of all available Pe and CCD times of light minimum were computed with 
the linear ephemeris given by Kreiner et al. (2001), 



Min.I = JX^eZ. )2428648.2837+ r.45448686 x E, 



(1) 



where 2428648.2837 is the time of conjunction, l d . 45448686 the constant ephemeris period 
P e , and E the cycle. The (O — C)i values are listed in the sixth column of Table 1 and are 
represented graphically against E in the upper panel of Figure 1. 

If the orbital period of a binary varies linearly in time, the (O-C) diagram should have 
a parabolic shape, as observed in other eclipsing binaries (e.g., Qian & Zhu 2006). As for 
V Puppis, a simple second-order polynomial does not give a satisfactory fit (the dashed line 
in the upper panel of Figure 1). Therefore, a sinusoidal term was added to a quadratic 
ephemeris to give a good fit to the (O — C)i curve. A least-squares solution leads to the 
following ephemeris, 



The sinusoidal term suggests an oscillation with a period of about T = 5.47 year and an 
amplitude of about A = 0. d 0163, which suggests a cyclic period oscillation with a period of 
5.47 years (the upper panel of Fig. 1). Though the (0 — C) 1 data do not cover a whole cycle, 
as shown in Figure 1, the minimum, the maximum, and the period of the cyclic variation 
are all well constrained, indicating that the period changes are reliable. The quadratic term 
in Eq. (2) reveals a continuous period increase at a rate of dP/dt = +2.10 x 10~ 7 days/year 
corresponding to a period increase of 1.8 s per century. 



Both components of V Puppis are B-type stars containing a convective core and a 
radiative envelope, which is different from the Sun. Hence, the period oscillation can not be 
explained as solar-type magnetic activity cycles (i.e. the Applegate mechanism) (Applegate 
1992). Therefore, the simplest explanation of the cyclic period oscillation is a wobble of 
the binary's barycentre due to the presence of a tertiary companion. The (O — C)2 values 
calculated from the quadratic ephemeris in Eq. (2) are shown in the seventh column of 
Table 1. The corresponding (O — 6% diagram is plotted against the epoch number E in the 
middle panel of Figure 1. By considering a general case with an elliptical orbit, the following 
equation was used to describe those (O — C)i residuals, 



Min.I = 2428648.2960(±0.0045) + l d .45448114(±0.00000091) x E 
+4.18(±0.44) x 10~ 10 x E 2 

+0.0163(±0.0016) sin[0.°2621(±0.0011) x E + 44.°3(±3.°9)]. 



(2) 



3. Discussion of the orbital period variation 



3.1. The presence of a massive black hole companion 



2 




(3) 
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The results were obtained with the same value of £1 (0.°2621) as that used in equation (2) 
(solid line in the middle panel of Figure 1). The orbital parameters of the tertiary companion 
were computed by using the formulae derived by Kopal (1959), 

a' 12 sin i' = cJa\ + bj, (4) 



, {b\ - a\)b 2 + 2a 1 a 2 b 1 

uu = arctan — — — -, (6) 

(a\ - b\)a 2 + 2a l b l b 2 

i , T a x b 2 - b x a 2 

t = t - — arctan — — , 7 

In a^a 2 + b\b 2 

where c is the speed of light, a' 12 , i', e', u>', and r' are the semi-major axis, the orbital 
inclination, the eccentricity, the longitude of the periastron from the ascending node, and 
the time of the periastron passage, respectively The solutions are displayed in Table 2. With 
the following well-known equations, 

An 2 

f(m) = x (a 12 sim') 3 (8) 

and 

f(m) = (M3Sin ^ 3 (9) 
J[ ] (Mi + M 2 + M 3 ) 2 ' [) 

the mass function, the masses and radii of the third body were calculated using the absolute 
parameters determined by Andersen et al. (1983). In the formula, Mi, M 2 , and M3 are 
the masses of the eclipsing pair and the third body, respectively, G is the gravitational 
constant, and T is the period of the (O-C) oscillation. For different values of the orbital 
inclination, the masses and orbital radii of the tertiary component are displayed in Table 
2. The corresponding relations between the mass and the radius of the third body and the 
orbital inclination are shown in Figure 2. 

The lowest mass of the tertiary companion is 10.4 M , i.e. larger than the mass of the 
secondary component. If the tertiary component is a normal star, we would see its spectral 
lines not changing with the orbital phase of the binary. However, this has not been reported 
(e.g., York et al. 1976; Koch et al. 1981; Andersen et al. 1983). Moreover, the third body 
is more massive than the secondary component (M 2 = 7.76 M ) in the central binary, it 
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should be very luminous and should contribute a large amount of third light to the total 
system unless it is not a normal main-sequence star. Therefore, we suspect that the tertiary 
component is a compact object and a probable black hole candidate. This situation resembles 
that in the triple system HR 2876 where a compact object (possibly a white dwarf star) with 
a mass excess of 1.0 M Q was found to be orbiting a B-type close binary star (B3.5V+B6V; 
P=15 days) (Burleigh & Barstow 1998; Vennes 2000). The main differences between the two 
systems are that the tertiary companion in the V Puppis is much more massive and that the 
central eclipsing binary in V Puppis system has a much more tight orbit and shows strong 
interaction between both components. 

A weak X-ray source (3U 0750-49) with 9.4 ± 2.3 count s -1 was detected by the Uhuru 
satellite (Giacconi et al. 1974). The counterpart of the X-ray source was a puzzling problem 
for astronomers since the 70s of last century. V Puppis was initially listed in the 3U catalogue 
as a possible candidate by Giacconi et al. (1974). However, another bright star HD 64740, 
close to the source (within the 90 percent confidence error box) was also considered as a 
possible candidate (Groote et al. 1978). By analyzing the X-ray data of V Puppis and HD 
64740 from the Copernicus satellite, Bahcall et al. (1975) later concluded that the weak 
X-ray source must almost certainly be closer to V Puppis than to HD 64740, but most 
probably neither V Puppis nor HD 64740 is the candidate. More recently, V Puppis was 
also discovered to be a X-ray source by ROSAT satellite with 0.08 count s" 1 (X-ray 50 A, 
250 eV) (e.g., Lampton et al. 1996). Our study indicates that the tertiary companion may 
correspond to the weak X-ray source. There are two mechanisms that can produce the X- 
ray emission in the V Puppis system. If the massive tertiary companion is a normal star, 
the X-ray emission could be caused by the colliding winds from the eclipsing pair and the 
tertiary companion, resembling those observed in other massive binary systems (e.g., Stevens 
1992; Pittard & Stevens 2002). If the third body is a massive black hole, the X-ray source 
is reproduced by accreting material from the massive binary via stellar wind as observed in 
X-ray emitting binaries (e.g., Prestwich et al. 2007). However, the former mechanism can 
be ruled out, because no spectral lines of the tertiary component were found. 

Another interesting feature of V Puppis is the existence of a significant amount of 
circumstellar matter, which was deduced from the optical spectra and ultraviolet observations 
obtained by the IUE spacecraft and the Copernicus satellite (e.g., York et al. 1976; Koch 
et al. 1981; Andersen et al. 1983). Numerous interstellar spectral lines, especially those 
of many types of heavy elements (e.g., Mg II, Si I, Si II, Si III, C I, O I, Fe II and Fe III) 
were found in the H II region around V Puppis, in which strong lines of Si III and Fe III are 
remarkable (e.g., York et al. 1976 and Koch et al. 1981). The formation of the circumstellar 
matter around V Puppis is an unsolved problem. Considering that V Puppis has undergone 
mass exchange, one may think that the circumstellar matter results from from the binary 



- 6- 



evolution. However, we may question why the other systems (e.g. AI Cru) that resemble V 
Puppis do not have the same significant amount of circumstellar matter (Bell et al. 1987). 
The most probable reason for the formation of the circumstellar matter around V Puppis 
is that it is formed by the supernova explosion of the progenitor of the tertiary black hole 
companion (e.g., Woosley & Weaver (1995) and Zhang et al. (2007)). This may produce the 
heavy elements in the circumstellar matter and results in the formation of the H II region 
around the eclipsing binary. 



3.2. The slow mass transfer of V Puppis 

V Puppis is a semi-detached binary system with the secondary component filling the 
critical Roche lobe (e.g., Andersen et al. 1983). The long-term period increase can be 
explained by the mass transfer from the secondary to the primary component. By using the 
equation, 

r N = 10 10 M 2 /L 2 , (10) 

where M 2 and L 2 are the masses and the luminosity of the less massive component, the 
nuclear time-scale of the secondary is calculated to be tn = 7.1 x 10 6 years. The value is 
close to the timescale of the period change r P = £ = 6.9 x 10 6 years, which reveals that V 
Puppis is now undergoing a slow mass-transfer evolutionary stage on the nuclear time-scale 
of the less massive component. This suggests that V Puppis is formed via a case A evolution 
and has passed through a rapid mass-transfer stage. The results are in agreement with the 
prediction from the theory of massive binary evolution (e.g., Sybesma 1986). By assuming 
a conservative mass transfer, the well-known equation 

P -1 1 

- = 3M(- -, (11) 

P V M 2 Mi v ; 

yields a rate of mass transfer dM/dt = 7.82 x 10~ 7 M Q /year. However, since the nuclear 
time-scale of the secondary is slightly longer than the time scale of period change, the slow 
mass transfer is insufficient to cause the observed period increase. This suggests that the 
stellar wind from the massive component of the binary should contribute to the period 
increase. The situation resembles that of AI Cru (Qian et al. 2008). 

By considering conservative mass and angular momentum, the calculation by Plavec 
(1968) suggests that, to produce a case A mass transfer of massive binaries, the initial 
orbital period should less than 1.8 days. One may ask a question what is the reason that 
causes the origin of the initially short-period detached system? It is possible that the massive 
tertiary companion has played an important role for the origin and evolution of inner binary 
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by removing angular momentum from the central system via Kozai oscillation (Kozai 1962) 
or a combination of Kozai cycle and tidal friction (e.g., Fabrycky & Tremaine 2007). This 
makes the inner eclipsing pair to have a lower angular momentum and a shorter initial orbital 
period. 

4. Conclusions 

The cyclic change of the orbital period, the presence of a weak X-ray source near V 
Puppis, and the significant amount of circumstellar matter around the binary star all support 
the conclusion that a massive black hole is orbiting the massive close binary V Puppis with 
a period of 5.47 years. Several massive stellar black hole candidates have been discovered 
recently (e.g., Greiner et al. (2001) and Orosz (2007)), but all candidates are found to 
be a component in a binary system. We conclude here that we have provided observational 
evidence for the first massive stellar black hole candidate as a tertiary companion of a massive 
close binary star. Theoretical investigations have shown that the formation of black hole 
binaries should experience a common-envelope phase with the spiral-in of the companion 
in the envelope causing the ejection of the envelope (e.g., Taam & Sandquist (2000) and 
Podsiadlowski et al. (2002)). The mean distance between the black hole candidate and 
the eclipsing pair is about 5.5 astronomical units (AU). During the supergiant phase, the 
radius of the precursor of the black hole companion in V Puppis may reach to several AUs. 
The survival of V Puppis in this evolutionary stage suggests that the tertiary black hole 
companion must be formed through an evolutionary process that is different from that of 
the black holes in X-ray binaries (e.g., Woosley & Weaver (1995) and Zhang et al. (2007)) 
and it could give some constraints to the evolution of single massive stars and the formation 
of isolated black holes. 

From the present investigation, we deduce that the progenitor of V Puppis is a triple 
system that is composed of a short-period detached B-type eclipsing binary and a much 
more massive tertiary companion. By drawing angular momentum from the inner binary 
through Kozai oscillation (Kozai 1962) or a combination of Kozai cycle and tidal friction (e.g., 
Fabrycky & Tremaine 2007), the massive tertiary companion may play an important role for 
the origin the central system, which causes the inner eclipsing binary to have a very short 
initial orbital period. In the triple system, the massive tertiary companion evolves faster, and 
finally into a massive black hole via a supernova explosion (e.g., Woosley & Weaver (1995) 
and Zhang et al. (2008)). At the same time, the original more massive component of the 
eclipsing pair evolves to fill its critical Roche lobe. Then, the binary was undergoing a rapid 
mass-transfer stage, and after the mass ratio was revered, the central binary finally reached 
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the present semi-detached configuration with a slow secondary-to-primary mass transfer on 
the nuclear time scale of the present less massive star. Supergiant stars have recently been 
found to be tertiary companions of OB-type close binaries in the two triple systems HD 
167971 and FR Set (e.g., Leitherer et al. (1987); Davidge & Forbes (1988), and Pigulski & 
Michalska (2007)). They resemble the progenitor of V Puppis with the tertiary supergiant 
companion evolving into a compact object. As for V Puppis, the existence of the black hole 
in the outer orbit poses strong constrain on the formation of the black hole, the mass loss in 
the supernova that lead to the formation of the black hole, and the magnitude of the kick 
received by the black hole, all make the object a much more interesting system for future 
study. 
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Table 1: All available photoelectric and CCD times of light minimum for V Puppis. 
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Table 2: Calculated orbital parameters of the tertiary companion in V Puppis. 
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Fig. 1. — A plot of the (O — C)i curve of V Puppis with respect to the linear ephemeris 
given by Kreiner et al. (2001). The solid line in the upper panel suggests a combination of a 
long-term period increase and a cyclic change, while the dashed line refers to the continuous 
increase of the orbital period. The middle panel shows the plot of the (O — C) 2 residuals 
of the eclipse timings based on the quadratic ephemeris in Eq. (1). The solid line in the 
middle panel is the theoretical orbit of the massive black hole candidate with an eccentricity 
of 0.45. Residuals after subtracting all effects of period change are displayed in the lower 
panel, and no variations can be traced there. 
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Fig. 2. — The relations between the mass M3 (M ) and the orbital radius d% (AU) of the 
tertiary component and its orbital inclination i' in the V Puppis system. The tertiary 
companion should be more massive than the secondary component and thus may be a black 
hole candidate. 



